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The freeze-out curve, which describes a vast amount of precise experi- 
mental data in heavy ion collisions, provides a relation between the colliding 
energy and the thermodynamical parameters of the fireball. The variance, 
skew and kurtosis of the event distribution of baryon number are studied at 
several energies of interest through Pade resummations of our Lattice QCD 
results. A smooth behaviour is predicted for three ratios of these quantities 
at current RHIC and future LHC energies. Any deviations from these at 
the RHIC energy scan would signal the presence of a nearby critical point. 
Our lattice results on the critical point do show such a behaviour. 
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1. Introduction 

Critical points in a phase diagram in the temperature-density plane are 
special for many reasons. Universality of critical indices, diverging corre- 
lation length are some of them. For common substances, such as water 
or carbon dioxide, the existence of critical point has been established ex- 
perimentally, with its location known rather precisely. For many of these, 
however, getting a theoretical, especially first principles based, computation 
of their locations turned out to be illusive still. I wish to describe to you 
how things are different for strongly interacting matter, which is naturally 
described by Quantum Chromo Dynamics (QCD) and has an inherently 
higher energy scale. Whether the QCD phase diagram has a critical point 
in its T-/j,b plane, where \xb is the baryonic chemical potential, is therefore 
a tougher question to address experimentally. Thanks to the impressive de- 
velopments on the experimental fronts, the Relativistic Heavy Ion Collider 
(RHIC) at BNL, New York, and the upcoming facilities such as FAIR at 
GSI, Darmstadt and NICA at Dubna, Russia can potentially search for it, 
as these proceedings reveal. 
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A variety of models have been successfully tested for hadronic interac- 
tions. These were was also the first set of tools used for getting a glimpse of 
the QCD phase diagram. For instance, using an effective chiral Nambu- Jana 
Lasinio type model, a phase diagram was obtained, which suggests pQ a crit- 
ical point to exist in a world with two light quarks and one heavier quark. 
One would clearly like to have an ab initio theoretical evidence for it. This 
turns out to be difficult as one usually has to deal with large coupling con- 
stants in the world of (low energy) hadronic interactions. Non-perturbative 
Lattice QCD, defined on a discrete space-time lattice, has proved itself to 
be the most reliable technique for extracting such information from QCD. 
The hadron spectrum has been computed and predictions of weak decay 
constants of heavy mesons have been made. Application of this approach 
to finite temperature QCD has yielded a slew of thermodynamics determi- 
nations, such as the pressure as a function of temperature. It is therefore 
natural to ask whether lattice QCD can help us in locating the QCD crit- 
ical point. The prospect of finding the critical point experimentally makes 
it exciting both as a check of theoretical predictions and as a competition 
for getting there first. Clearly, in view of the complexity of this task, one 
could turn to lattice QCD again to see if it can provide any hints for the 
experimental search program. In this talk, I summarise the results of our 
group in TIFR, Mumbai in these directions. 

Due to the well-known fermion doubling problem, one has to make a 
compromise in choosing the quark type for any computation. Mostly stag- 
gered quarks are used in the lattice simulations at finite temperature and 
density. These have an exact chiral symmetry which provides an order pa- 
rameter for the entire T-fiB plane but unfortunately the flavour and spin 
symmetry is broken for them on lattice. A representation of precisely 2 
(or 2 +1) quark flavours may thus be problematic, more so on the coarse 
lattices one is constrained to employ. The existence of the critical point, on 
the other hand, is expected to depend crucially on the number of flavours. 
Although computationally much more expensive, domain wall or overlap 
fermions are better in this regard, as they do have the correct symmetries 
for any lattice spacing at zero temperature and density . Introduction of 
chemical potential, fj,, for these, however, is not straight-forward due to their 
non- locality. Bloch and Wettig [2] proposed a way to do this. Unfortunately, 
it turns out [3] that their prescription breaks chiral symmetry . Further- 
more, the chiral anomaly for it depends on [i unlike in continuum QCD [4 J . 
While both these deficiencies go away in the continuum limit, the practical 
problem of the lack of an order parameter on finite lattices remains. What 
is therefore desperately needed is a formalism with continuum-like (flavour 
and spin) symmetries for quarks at nonzero [i and T with a well-defined 
order parameter on lattice. 
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Finite density simulations needed for locating a critical point suffer from 
another well known problem. This one is inherited from the continuum the- 
ory itself: the fermion sign problem. It is a major stumbling block in ex- 
tending the lattice techniques to the entire T-\xb plane. Several approaches 
have been proposed in the past decade to deal with it. Let me provide 
a partial list: 1) Two parameter re-weighting [5], 2) Imaginary chemical 
potential [6\, 3) Taylor expansion [7J, 4) Canonical ensemble method [3J, 
and 5) Complex Langevin approach [D|. We employ the Taylor expansion 
approach, developed independently by the TIFR group [7| and the Bielefeld 
group [7], to obtain the results discussed in the next section. 



2. Lattice Results 

Our results were obtained by simulating full QCD with two flavours of 
staggered fermions of mass m/T c = 0.1 on Nt x N^ lattices, with Nt = 4 
and N s = 8, 10, 12, 16, 24 [10J and a finer N t = 6 with N s = 12, 18, 
24 [TT]. Earlier work by the MILC collaboration for the smaller Nt lattice 
had determined m p /T c = 5.4 ± 0.2 and m n /m p = 0.31 ± 0.01, leading to a 
Goldstone pion of 230 MeV in our case. For the finer lattice, we determined 
/3 C as well. We covered a temperature range of 0.89 < T/T c < 1.92 by 
suitably choosing the range of couplings on both lattices. Our measurements 
were made typically on 50-200 independent configurations, separated by the 
respective autocorrelation lengths. 

From the canonical definitions of number densities and susceptibilities 
Xn u ,n d -, the QCD pressure P can be seen to have the following expansion in 
p. 

AP _ P(fx,T) P(0,T) v J_ ff^y- J_ ffM\ nd , 91 . 

T 4 - T 4 T 4 A.Xn u ,n d n j y T J n d \ \T ) 1 ' 

where the indices n u and denote the number of derivatives of the par- 
tition function with respect to the corresponding chemical potentials. We 
construct the series for baryonic susceptibility from this expansion, and look 
for its radius of convergence as the estimate of the nearest critical point. 
Successive estimates for the radius of convergence were obtained by using 

the ratio method [fn+l,n+3 

= Jn(n + 1)Xb +1) /Xb +3)t2 ] and the n ^ root 

/ (2) (n+2) \l/ n 

method [r n = (n\x B /Xb T n \ ]. We used terms up to 8th order in \x 
for doing so. A key point to note is that all coefficients of the series must 
be positive for the critical point to be at real [a, and thus physical. We 
thus first look for this condition to be satisfied and then look for agreement 
between the two definitions above as well as their n-independence to locate 
the critical point. The detailed expressions for all the terms can be found in 
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|10j where the method to evaluate them is also explained. We use stochastic 
estimators. For terms up to the 8 th order one needs 20 inversions of (D + m) 
on ~ 500 vectors for a single measurement on a given gauge configuration. 
This makes the computation expensive. Nevertheless, extension to 10th 
& even 12th order may be possible, especially using new emerging ideas [4] 
which may save up to 60 % computer time in the measurements. Adding the 
chemical potential simply as (j,N, as proposed in [I], one has the advantage 
that most derivatives of the quark matrix with \i are zero except the first. 
This reduces the number of terms in each nonlinear susceptibility (NLS) 
appearing in the equation above a lot. One can show that it still leads to 
essentially the same results [12]. 




Fig. 1. Estimates of radii of convergence as a function of the order n and for our 
two methods at T/Tc = 0.99 and the critical point T E /T c =0.94. The open (filled) 
points display the results for the ratio (root) method. From Ref. |llj . 



Fig. [T] shows our results for both the ratios defined above on JV; = 6 
lattice at two different temperatures, T/T c = 0.99 and 0.94. All the suscep- 
tibilities are positive at both but the ratios fluctuate for the former while 
seeming to be constant for the latter. WepT] thus found the coordinates 
of the endpoint (E) — the critical point — to be T E /T c = 0.94 ± 0.01, and 
fj, E >/T E = 1.8 ± 0.1 for the finer lattice. Our earlier coarser [TU] lattice 
result was ^/T E = 1.3 ± 0.3 for similar volumes with an infinite volume 
result leading to H%/T E = 1.1 ± 0.1. In view of this, it may be safer to 
quote the critical point to be at /j,b/T ~ 1 — 2. 

As a cross check on the location of the n E /T E ', one can use the series 
with the numerically determined coefficients directly to construct xb for 
nonzero /j,. As one sees in Fig. [2j it leads to smooth curves with no signs 
of criticality, whereas employing the Pade approximants for the series to 
estimate the radius of convergence does lead to a consistent window with 
our above estimate. 
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Fig. 2. The baryonic susceptibility as a function of the baryonic chemical potential 
using simple polynomial with all computed orders (left) and Pade approximants 
(right). From Ref. [TT] . 

3. Searching Experimentally 

In order to make predictions for the heavy ion experiments, and design 
further any search criteria for the critical point, one needs to know the 
temperature T and the chemical potential /i^ as a function of the colliding 
energy. Recently, we proposed [13] to use the freeze-out curve as a a way to 
exploit the information hidden in the nonlinear susceptibilities discussed in 
the previous section by evaluating lattice QCD predictions along it. Recall 
the well known observation that the huge amounts of precise data on hadron 
abundances in a variety of relativistic heavy ion experiments are well de- 
scribed [13] by statistical thermodynamical models. Indeed this description 
has been found to work over the entire range of experiments made so far. 
It thus leads to a mapping of the T and \xb parameters of the model to 
the collision (center of mass) energy \fs. Assuming these to reflect the true 
thermodynamic variables of the system and plotting these results in the T- 
\xb plane, one has the freeze-out curve. We treat this curve, solely based on 
precise experimental data on hadron yields, as a means to provide the (T, /x) 
accessible in heavy-ion experiments as a function of the colliding energy ^/s. 
We use then these relevant freeze-out parameters T and in our lattice 
computations directly to make predictions as a function of colliding energy. 

One expects the QCD critical point to have a critical region whose 
size depends on the size of the fireball as well as the critical exponents. 
The freeze-out curve may pass through it or may miss it. If the condi- 
tions are thus right, it may pass close enough to the critical point such 
that a study of fluctuations along it will detect its presence. We define 
mi = Tx^Ht^b)/^^,^), m 3 = T X ^(T, n B )/x®(T, »b), and m 2 = 
minis. Note that the variance, skew and kurtosis of the event distribu- 
tion of baryon number measure the various x' s above. Furthermore, the 
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spatial volume cancels out in these ratios, making them also suitable for 
experiments who can use their favourite proxy for it. Usually, number of 
participants is preferred for that role in the analysis of heavy ion data. 

Defining z = hb/T, and denoting by ry the estimate for radius of con- 
vergence using Xi, Xj as before, one has 
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Similar series expressions can be written [13] down for 1712 and 7713. We 
resum these by the Pade method described above to construct them, since 
they seem to capture the critical behaviour well in Fig. [2j 
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Fig. 3. Results for m, on the N t =6 lattice as a function of the colliding energy. A 
power law fit to the high energy data and a hadron resonance gas fit is also shown. 
FromRef. [13]. 



Figure [3] shows our [13] results for these ratios on our finer lattice; the 
results for coarser lattice are similar and can be found in Ref. [13]. Also 
shown in these figures are hadron resonance gas results which do not as- 
sume a critical point. One observes a smooth and monotonic behaviour for 
large -y/i which is well reproduced by the hadron resonance gas. Note that 
even in this smooth region, any experimental comparison is exciting since 
it is a direct non-perturbative test of QCD in hot and dense environment. 
Earlier lattice predictions, such as the equation of state or the transition 
temperature can be compared with experiments only indirectly. Remark- 
ably, a non-monotonic behaviour is visible in Fig. Oat our estimated critical 
point in all mj, suggesting that it would be accessible to the low energy scan 
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of RHIC at BNL. Indeed, even if one were to be cautious in trusting the 
numerical precision of our results, what should be clear that such a non- 
monotonic behaviour seen at any other nearby location would still signal 
the presence of QCD critical point. 

In order to confront these results on the baryon number fluctuations 
with data, one needs to address the issue of neutral baryons — neutrons — 
which are not easy to detect and are thus missed. It turns out that proton 
number fluctuations suffice |15j . Since the diverging correlation length at 
our critical point is linked to the a mode which cannot mix with any isospin 
modes, and the isospin susceptibility xi must be regular there. Assuming 
protons, neutrons, pions to dominate, Ref [15] showed xb to be dominated 
by proton number fluctuations only. The STAR collaboration has recently 
exploited this idea and constructed the ratios mi and 777,2 from net proton 
distributions [16J. 
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Fig. 4. Comparison of the STAR collaboration results with our mi results as a 
function of number of participants. From Ref. [16] . 

Figure d] shows their results for mi against our lattice data. Remarkably, 
one observes a good agreement with our lattice results. Hopefully, future 
results from the RHIC energy scan will locate the critical point this way. 
Similar agreement is also seen for 7772. 



4. Summary 

The QCD phase diagram in the T-/j,b plane has begun to emerge using 
first principles lattice approach. Our finer lattice results for Nt = 6 are first 
to begin the crawling towards continuum limit, suggesting a critical point 
at Hb/T ~ 1-2. Ratios of nonlinear susceptibilities appear to be smooth 
on the freeze-out curve at large colliding energy. Critical point leads to 
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structures in these mj, which may be accessible in future. STAR results on 
proton number fluctuations appear to agree with our lattice QCD predic- 
tions, making this a unique direct non-perturbative test of hot and dense 
QCD in experimental tests. So far no signs of a critical point have emerged 
in the CERN and RHIC experimental results. An interesting question is 
whether the RHIC energy scan deliver it for us or whether we will need to 
wait for FAIR at GSI to be operational. At any rate, exciting future awaits 
these programs and all of us interested in it. 
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